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Abstract 



We demonstrate that pure space-like axial gauge quantizations of gauge fields 
can be constructed in ways which are free from infrared divergences. To do so we 
must extend the Hamiltonian formalism to include residual gauge fields. We construct 
an operator solution and an extended Hamiltonian of the pure space-like axial gauge 
Schwinger model. We begin by constructing an axial gauge formulation in auxiliary 
coordinates: x M = (x + ,x~), where x + = x° sin 9 + x 1 cos 9, x~ = x° cos 9 — x 1 sin 9 
and we take A_ = ^4° cos 9 + A 1 sin 9 = as the gauge fixing condition. In the region 
0<# < f , we can take x~ as the evolution parameter and construct a traditional canon- 
ical formulation of the temporal gauge Schwinger model in which residual gauge fields 
dependent only on x + are static canonical variables. Then we extrapolate the temporal 
gauge operator solution into the axial region, | < 6 < |, where x + is taken as the evo- 
lution parameter. In the axial region we find that we have to change representations of 
the residual gauge fields from one realizing the PV prescription to one realizing the ML 
prescription in order for the infrared divergences resulting from (<9-) _1 to be canceled 
by corresponding ones resulting from the inverse of the hyperbolic Laplace operator. 
We overcome the difficulty of constructing the Hamiltonian for the residual gauge fields 
by employing McCartor and Robertoson's method, which gives us a term integrated 
over x~ = constant. Finally, by taking the limit 9 — > | — we obtain an operator 
solution and the Hamiltonian of the axial gauge (Coulomb gauge ) Schwinger model 
in ordinary coordinates. That solution includes auxiliary fields and the representation 
space is of indefinite metric, providing further evidence that "physical" gauges are no 
more physical than "unphysical" gauges. 
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1. Introduction 



Axial gauges, n^A^ = 0, specified by a constant vector n M , have been used recently in 
spite of their lack of manifest Lorentz covariance. One reason is that the Faddeev-Popov 
ghosts decouple from the theory in the axial gauge formulations. 1 -* The case of n 2 = 0, the 
light-cone gauge, has been extensively used in light-front field theory (LFFT) in attempts 
to find nonperturbative solutions of QCD. In LFFT one usually makes the change of vari- 
ables, xf = ^-^p, xf = z ° v ^ 3 , 2 \ and specifies quantization conditions on the hyperplane 
x~l = 0. Advantages are that the vacuum state contains only of particles with nonnega- 
tive longitudinal momentum and that there exist relativistic bound-state equations of the 
Schrodinger-type. (For a good overview of LFFT see Ref. 3).) 

It has been known for some time that axial gauge formulations are not ghost free, contrary 
to what was originally expected and is still sometimes claimed. It was first pointed out by 
Nakanishi 4 ) that there exists an intrinsic difficulty in the axial gauge formulations so that 
an indefinite metric is indispensable even in QED. It was also noticed that in order to bring 
perturbative calculations done in the light-cone gauge into agreement with calculations done 
in covariant gauges, spurious singularities of the free gauge field propagator have to be 
regularized not as principal values (PV), but according to the Mandelstam-Leibbrandt (ML) 
prescription 5 ^. Shortly afterwards, Bassetto et al. 6 ^ obtained the ML form of the propagator 
in a canonical formalism. They quantized at equal time in light-cone gauge, which required 
the introduction of a Lagrange multiplier field and its conjugate. These may be viewed as 
residual gauge degrees of freedom. Furthermore, Morara and Soldati 7 ) found just recently 
that the same is true in the light-cone temporal gauge formulation in which x~t is again taken 
as the evolution parameter but A "J^ 3 = is taken as the gauge fixing condition. McCartor 
and Robertson 8 ^ showed that in the light-cone axial gauge formulation, where A °^' 4 = 
is taken as the gauge fixing condition, the translational generator P + consists of physical 
degrees of freedom integrated over the hyperplane xf = constant and residual gauge degrees 
of freedom integrated over the hyperplane xj = constant. 

In McCartor and Robertson's work the residual degrees of freedom were viewed as in- 
tegration constants necessary to completely specify the solution to the constraint equations 
which relate constrained degrees of freedom to the independent degrees of freedom. That 
is the attitude we will adopt in the present paper. Such equations often occur in quantum 
field theory and little attention has been given to the question of what boundary conditions 
are appropriate to completely specify their solution. In the present paper we will attend to 
that question. An example of the type of equation we have in mind is the equation for A + 
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in the light-cone gauge 

a 2 A- + <9_<m i = j_. (i-i) 

This equation admits a solution of the form 

C(x + , ar 1 ) + ar£(a;+ ar 1 ) + 

where F(x) is any solution to the equation. McCartor and Robertson showed that if F(x) 
is taken to be the usual solution, which results from the replacement <9 M — > ik^ in the mode 
expansion, then neither C(x + ,x ± ) or B(x + ,x ± ) are zero in the case of QED. On the other 
hand, in an an earlier work 9 -* , we showed that for the Schwinger model in the light-cone 
gauge, B(x + ,x ± ) is zero while C(x + ,x ± ) is not. It will be seen below that it is this last 
case that is peculiar: for all the solutions we shall give, both B(x + ,x ± ) and C(x + ,x ± ) 
are nonzero except for the special case of the light-cone gauge where B(x + , x^) happens to 
vanish. 

Because the axial gauges can be viewed as continuous deformations of the light-cone 
gauge, extensions of the ML prescription outside the light-cone gauge formulation have also 
been studied. Lazzizzera 10 ** , and Landshoff and Nieuwenhuizen 11 - 1 , constructed canonical for- 
mulations for the non-pure space-like case (n°^0, n 2 < 0) in ordinary coordinates. However, 
in spite of many attempts, no one has succeeded in constructing a consistent, pure space-like 
axial gauge (n° = 0, n 2 < 0) formulation. One problem is that when the residual gauge 
fields are introduced as integration constants that occur when one integrates the axial gauge 
constraints, one cannot obtain their quantization conditions from the Dirac quantization 
procedure 12 - 1 . Another problem is that the residual gauge fields are independent of x 3 in the 
gauge A3 = so that the Hamiltonian cannot be obtained by integrating densities made of 
those residual gauge fields over the three dimensional hyperplane x° = constant. 

This motivates us to study the problem of introducing the residual gauge fields into 
pure space-like axial gauge formulations. To construct the Hamiltonians and to determine 
quantization conditions, we construct the pure space-like gauge as a continuous deformation 
of the light-cone gauge. In a previous, preliminary work 13 ) we have constructed an axial 
gauge formulation of noninteracting abelian gauge fields in the auxiliary coordinates x M = 
(x + , x~ , x 1 , x 2 ), where 

x + = x sin# + x 3 cos#, x~ = x°cos6 — x 3 sin8 (1-2) 

and A_ = A°cos9 + A 3 sin8 = is taken as the gauge fixing condition. The same framework 
was used previously by others to analyze two-dimensional models. 14 ^ In the region 0<9 < | 
we take x~ as the evolution parameter and construct the canonical temporal gauge formu- 
lation; in that case the residual gauge fields are static canonical variables. By continuation, 
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we can obtain operators satisfying the field equations in the axial region ^ < 9 < |, where 
x + is taken as the evolution parameter. In the axial region we cannot use the traditional 
way of constructing the Hamitonian for the residual gauge fields; we obtain it by integrat- 
ing the divergence equation of the energy-momentum tensor over a suitable closed surface. 
Because the residual gauge fields do not depend on x~ , their boundary surface (x~ — > ±00) 
contributions have to be kept. As a consequence, we obtain the Hamiltonian which includes 
a part from integrating a density involving the residual gauge fields over x~ = constant. 
Then, by taking the limit 9 — > | — we obtain the Hamiltonian in the axial gauge A 3 = 
formulation in ordinary coordinates. 

In this paper we proceed to interacting gauge fields. As a first step in finding ways 
to introduce the residual gauge fields into interacting axial gauge theories, we consider the 
exactly solvable Schwinger model. Previous attempts to construct the axial gauge (Ai = 0) 
operator solution to the Schwinger model have used a representation space isomorphic to 
(copies of) that of a positive metric free massive scalar field E. In other words, they have 
assumed that {A\ = 0) is a physical gauge in that the entire representation space is physical. 
All such attempts encounter severe infrared difficulties. 15 -'. We expect that the difficulties 
can be overcome by introducing residual gauge fields depending only on x°. To study that 
possibility, we consider constructing an operator solution of the pure space-like axial gauge 
Schwinger model as a continuation from the light-cone gauge. We begin by constructing an 
operator solution in the region 0<9 < |. This is straightforward since we can use the axial 
anomaly and the temporal gauge canonical quantization conditions as guiding principles in 
constructing the operator solution. We continue the operator solution into the axial region 
I < 9 < |. In the axial region the operator (<9_) _1 becomes singular and gives rise to 
infrared divergences. Furthermore, in the axial region the Laplace operator, m 2 — n_<9 + 2 
(where m 2 = — and n_ = cos 29), becomes hyperbolic so that divergences result from its 
inversion. We find that the former divergences are canceled by the latter ones if we choose 
the representation of the residual gauge fields in such a way as to realize the ML prescription. 
As a consequence we recover well-defined fermion field operators. Bilinear products of these 
give rise to bosonized expressions for the vector current and the energy-momentum tensor 
through the equal x + -time point-splitting procedure. We apply McCartor and Robertson's 
method 8 ^ and obtain the Hamiltonian for the residual gauge fields. It consists of a part 
from the physical operators integrated over x + = constant and a part from the residual 
gauge operators integrated over x~ = constant. This is the way we extend the traditional 
Hamiltonian formulation to obtain the Hamiltonian for the residual gauge fields. Now it is 
straightforward to obtain an operator solution and the Hamiltonian in the gauge A 1 — 
(here, the Coulomb gauge) in the ordinary coordinates. All one has to do is to take the 
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limit 9 — > f — 0. While that limit exists, we cannot simply set 9 equal to |. For that 
reason one might say that we have not really quantized the Schwinger model at equal time 
in the Coulomb gauge. On the other hand one might regard 9 as a regulating parameter. 
In that view we see that, as was the case in the light-cone gauge 9,16 ) , proper regulation of 
the theory necessarily involves information off the initial value surface. In any event, the 
solution contains the auxiliary fields and the representation space is of indefinite metric. This 
provides further indication that "physical" gauges are no more physical than "unphysical" 
gauges in that the representation space contains unphysical states and is of indefinite metric 
in all cases. 

The paper is organized as follows: In § 2, we use the axial anomaly and the temporal 
gauge quantization conditions as guiding principles to construct the temporal gauge operator 
solution of the Schwinger model in the auxiliary coordinates. We also point out that if we 
choose the representation of the residual gauge fields so as to realize the PV prescription, 
we can construct fermion field operators in such a way that vacuum expectation values of 
them vanish. In § 3, we continue the solution into the axial region and show in detail that 
the infrared divergences inherent to the axial gauge quantizations are canceled if we choose 
the representation of the residual gauge fields which realizes the ML prescription. We also 
discuss properties of the axial gauge solution. Section 4 is devoted to concluding remarks. 

We use the following conventions: 

g =cos26>, g_ + = g + _ = sin 29, g ++ = —cos29, 

g— = cos 29, g~ + = g + ~ = sin 29, g ++ = - cos 29, 
= (n + ,n~) = (0, 1), n M = (n+,n_) = (sin26>, cos26>), 
7° = <7i, 7 1 = ia 2 , 7 5 = -^3, 



7 + = 7 sin 9 + 7 1 cos 9 



+ 

—e(n-)y/l — n + 



a 1 • a ( £ ( n -)v /1 - ™4 

7 = 7 cos 9 — 7 sin 9 — , 

§2. Temporal gauge formulation in the auxiliary coordinates 

In this section we confine ourselves to the region 0<6> < | and choose x~ as the evolution 
parameter. The temporal gauge Schwinger model in the auxiliary coordinates is defined by 
the Lagrangian 

L = --F^ - B{n- A) + ifrf^ + ieA^ (2-1) 
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where B is the Lagrange multiplier field, that is, the Nakanishi-Lautrup field in noncovariant 
formulations. 17 ^ From the Lagrangian we derive the field equations 

d^ u = n u B + r, J v = efrfl> ( 2 - 2 ) 

i-f{d„ + ieA v )tl) = Q, (2-3) 

and the gauge fixing condition 

A_ = 0. (2-4) 

The field equation of B, 

d-B = 0, (2-5) 

is obtained by operating on (|2-2|) with d v . Canonically conjugate momenta are defined by 



8L , 5L 5L 5L - _ . , 

*~ = S^Al = °' * = 5dJ£ = F ^* B= 8W = ">** = 8W =x ' h ~- (2 ' 6) 

Thus we impose the following equal x~-time canonical quantization conditions: 

[A + (x),A+(y)} = [7r + (x),7r+(y)] = 0, [A + (x), ^{y)\ = i8(x+-y + ), (2-7) 

[A+(x)Mv)] = [k + (x)Mv)] = 0, (2-8) 

{M*),My)} = °» {Mz),r P (v)} = ((M" 1 )^-^). (2-9) 

We begin constructing an operator solution by imposing the chiral anomaly 

<9_ J+ - <9+J_ = -m 2 F_ + . (2-10) 
The current operator given by Q2-2p is 

J_ = g _ + J+ + g __J- = -d + F_ + -n_B, J + = d-F_ + -n + B. (2-11) 
Substituting these expressions into (|2-10|) yields 



(d^ + m 2 )F_+ = -n_d+B. (2-12) 
Because cLB = 0, F_ + is solved to be 



F_ + = cLA+ = -mS = 7^ d+B. (2-13) 



71- 

Here we have identified the homogeneous solution with the massive free field S, since F_ + 
is gauge invariant. The normalization is determined in such a way that it agrees with the 
Landau gauge solution 18 ). E is given by 

E{x) = -L T *^{ a (p + ) e -*- + a*(p + )e**}, (2-14) 

A/47T J-oo \/» _ 



where p = yp+ + rn 2 with m 2 = n_m 2 and 

[a(p+),a(q+)] = 0, [a(p+), a*(q + )\ = 5(p + - q+). (2-15) 
By integrating ( [2 -131) with respect to x~ we obtain the general solution for A + as 



TXXi ~ TXi CC 

A + = — — E tt^B + integration constant. (2-16) 

cL m 2 — n_o + 

We determine the integration constant in the following way: From the form of the d'Alembertian 
in the auxiliary coordinates, 

= n_d_ 2 + 2n + n„d+d- - n_d + 2 (2-17) 
we rewrite the first term of (|2-16 ) as 



cL m 2 — n^o + 

In the range of 9 we are considering, we see that the operator (<9_) _1 does not give rise to any 
infrared divergences. Now the first term in ( |2-16|) has the property [^-E(x + ), ^ L E(y + )} ^ 



at equal x -time because the E satisfies the following equal x -time commutation relations 
[E(x),E(y)} = [d-Z(x),d-£{y)} = 0, [E(x), drE{y)\ = i5(x + - y + ). (2-19) 



Therefore, to preserve the canonical commutation relations fl2-7|) , we must cancel the nonva- 
nishing terms by contributions from another field C, the integration constant, which depends 
only on x + . We find that if B and C satisfy the commutation relations 

[B{x),B(y)\ = [C{x),C{y)\ = 0, [B{x),C{y)\ = i{m 2 - n_d + 2 )5(x + - y+), 

[B{x),E(y)) = [C{x),E{y)} = 0, (2-20) 

then 

A + = + 2 9+ *2 { C - n -^ B ) + ( 2 n+m2 ,2, 2 B (2-21) 
cL m 2 — n_o + v J [mr — n_c* + j 2 

satisfies the quantization conditions in (|2-7| ). In this way we can determine the integration 
constant: by taking account of the canonical quantization conditions in the temporal gauge 
formulation. 

Now that we have the A+, we proceed to solve the fermion field equation Q2-3|) . A formal 
solution is given by 

V>a(z) = , Z ° exp[-iy/^A a (x)], (2-22) 
V(T°7")« Q 
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where Z a is normalization constant and 

AJx) = X(x) + (-l) a X(x) (2-23) 

with 

x = jT^ + 2 m *2 ( c - n - x ' B + 2 n+n " 2 9 + b] , 

cL m — n_o + \ m z — nJo+ ) 

\ = E + — j^d-'B. (2-24) 

mr — n_a + 

Here the operator 8+ 1 is defined by 

j poo 

d?B(x) = - / dy + e(x + - y + )B(y), (2-25) 

Z J —oo 

which imposes, in effect, the principal value regularization. To further examine the properties 
of the fermion field operators, such as their vacuum expectation values, we need explicit 
representations of the B and C. Note that the canonical commutation relations cannot 
uniquely determine the representations of the canonical variables, as was pointed out by 
Haller 19 -'. One possibility is to choose the following zero norm and positive metric expressions 
for B and C: 

B(x) = -= / dk + J\k^\{(a (k + ) - ai (k + ))e-' lk + x+ + h.c.}, 

V47T J -oo v 

2 m , 2 C(x) = ^=T ^={ai{k + )e- lk+x+ - al(k + y k +* + }, (2-26) 
mr — n-0+ -JAtt J-oo J\h . I 



where 



[a (k + ),a* (q + )} = -5(k + -q + ), [oi(fc+), oJ(g + )] = 5(k+ - q+), (2-27) 



and all other commutators are zero. It can be shown that B and C in Q2-26 ) realize the 



PV prescription. With this choice, the spurion operators are composed out of ultraviolet 
parts of the E as well as infrared parts of C and d^B according to the method given in 
Ref.18). This in turn gives rise to vanishing vacuum expectation values for the fermion 
field operators. In the Landau gauge solution 18 ) the fermion fields have vanishing vacuum 
expectation value but in the light-cone gauge solution the vacuum expectation value of ip + 
diverges. The inclusion of the ultraviolet parts of E in the spurions causes them to become 
x dependent and the solution has some different characteristics from the previous solutions. 
We are not completely certain whether the solution using the choice ( ^-2tj| ) is completely 
consistent but we do know that when we move to the axial gauge choice [6 > ?) we cannot 
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use the representation (|2-26|) : if we do there are uncancelled infrared divergences which 



destroy the solution. We will not further discuss the choice ( |2-26| ) in the present paper but 
will concentrate on another choice, one which realizes the ML prescription. To do that we 
must limit the £; + -integration region to be (0, oo); we showed in Ref.13) that by choosing a 
suitable vacuum it is always possible to do so. We show in the next section that with the 
new choice, the infrared divergences in the axial gauge cancel. 

To implement the new(ML) representation we choose the following representation for B 
and C: 



B(x) = -= / dk + Jkl{B(k + )e- ik+x+ + B*(k + )e ik+x }, 

V27T JO V 

' T C(x) = / ^L{C(k + )e~ ik +* + - C*(k + y k +* + }, (2-28) 

\/2tc Jo \ Au. 



m 2 — n_d+ y/2n Jo a/^+ 

where 

[B(k + ),C*(q + )] = [C(k + ),B*(q + )} = -S(k+ - q+), (2-29) 

and all other commutators are zero. Note that the C becomes a zero norm field as a 
consequence of altering the representations. Note also that the vacuum expectation value of 
ip2 diverges as it does in the light-cone gauge solution 9 ). 
We define the physical subspace, V, by 

V = { |phys) | £?(£; + )|phys> = }. (2-30) 

To subtract the infrared parts of the massless fields into the spurion operators a a , we define 
the infrared parts, of A a to be 

i f 1 ^ die 

4 0) = 4= / -7={C(k + ) - C*(k + ) + (~lT(B(k + ) - B*(k + ))} (2-31) 
\ 2tv Jo \/k-L 



/2n Jo \fk 

where k is a small positive constant. Then the fermion field operators are defined as 

il>Jx) = a exp[-iy/ZA(-\x)]a a exp[-iVZA£>(x)] (2-32) 

V / (7°7 _ )aa 

where A^) and A$ are the creation and annihilation operator parts of A ar =A a — A^ and 

Q 



<y a = exp 



-Wvr 



AS> - (-1, 2m 



(2-33) 



Here, Q = — /f^ dx + B(x); note that Q in a a constitutes a Klein transformation. 

We enumerate properties of the ip a to show that the ipa constitute the operator solution 
of the temporal gauge Schwinger model in the auxiliary coordinates. 
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(1) The Dirac equation is satisfied: 



i^(d^ + ieAJifi = 



(2-34) 



To get this we use the fact that [A + (x),ip a (x)] = 0, where A + is the annihilation operator 
part of A + . 

(2) The canonical commutation relations ( |2-8| ) and anticommutation relations (|2-9|) are 
satisfied. 

(3) By applying the gauge invariant point-splitting procedure to etp^ip, we obtain the 
vector current J M = me^ v d v \ where e ++ = e = 0, e + ~ = —e h = 1. The charge 
operator, Q, is given by 



where the derivative terms integrate to zero. 

(4) The conserved chiral current J^ 5 = me^ v d v X is similarly obtained from e-07 At 7 5, ?/>. 
The chiral charge operator, Qs, is given by 



where the derivative terms, except d + C, integrate to zero. The term d + C does not integrate 
to zero because ip contains the singular operator d+ l B. 

(5) Applying the gauge invariant point-splitting procedure to the fermi products in the 
symmetric energy-momentum tensor and subtracting a diverging c-number, which we denote 
by R, we get 




(2-35) 




(2-36) 



6 



+ 



-R(^-d + iJ - a^ 7 » - A + J~ - -(A + , B) + 

hd + \,d-\) + -hA + ,B) + , 



(2-37) 



e 



+ 



+ 




d + ^ + ^)-A + J + + -{F_ + ) 2 



+ 



(2-38) 



e 



--R^+d^ - d + ^ + ^) + A + J+ + i(F_ + ) 2 

, i 
-{(d + A) 2 + (^A) 2 } + -(F_ + ) 2 , 



(2-39) 
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2n_ 2n_ 

+ — J_A+ + ^± J + A + = -(<9_A, <9 + A) + . (2-40) 

n_ n_ 2 

(6) Translational generators consist of those of the constituent fields: 



•X- 



P_ = / dx + : OS : 

J — oo 



/°° I ~ ~ ~ 7? 

cfe+- : {n_(d + E) 2 + n_(d_£) 2 + m\Ef + £? = - 2 B} :, 
-oo 2 m — n_a. 



/CO /"CO _ _ j 
cta+ : - := / dx + : {d+Edr E - B— ^d + C} : . (2-41) 
-oo J-oo ' m 2 — n_o + 

§3. Cancelation of Infrared divergences resulting from dz 1 

In this section we continue 9 into the region j < 9 < | and take x + as the evolution 
parameter. In accordance with the change of the evolution parameter we change the Fourier 
expansion of E so as to be convenient for calculating equal x + -time commutation relations. 



By making use of = — -rp we change integration variables from p + to p_ and rewrite E 
as follows 



E(x) = -L f° ^{ a (p_)e-*- + a*(p_)e*-*}, (3-1) 



where p + = \JP- + rn 2 with m 2 = |n_|m 2 and 

[a(p_),a(<?-)] = 0, [a(p_),o*(g_)] = <J(p_ - <?_). (3-2) 
17 satisfies the following equal x + -time commutation relations: 

[E(x), E(y)} = [d + E(x),d + E(y)] = 0, [E(x),d+E(y)} = iS(x' - y~). (3-3) 

Furthermore, in accordance with the change of the evolution parameter, the conjugate mo- 
mentum of ifi changes to 

Therefore, we have to modify the normalization of ip in ( |2-32j ) accordingly: 

= I, Z f \, eM-iV^A^(x))a a exp[-iV^AW(x)], (3-5) 



where we have put (a) on ip and on the normalization constant to denote that the normal- 
ization is altered. 
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We begin by pointing out that because n_ = cos 29 < in the region | < 6 < |, 
the Laplace operator, m 2 — n_<9 + 2 , which operates on the residual gauge fields, becomes 
hyperbolic so that its inverse gives rise to divergences. Therefore, we regularize them by 
PV regularization. We then recover well-defined fermion field operators because linear in- 
frared divergences resulting from are canceled by those resulting from the square of the 
inverse of the hyperbolic Laplace operator. To show that this really happens, we rewrite 
ij}^* (x)i/j^ (y) into normal ordered form. It is tedious but straightforward to obtain 

V4 a) Wi a) (y) = \ a J exp[M aa (a; - y)]exp[-i^(A^(x) - A^\y))} 

{TT)aa 



x exv[-i^(AW(x) - AM(y))}, (3-6) 



where 




2 + 2(-l 




2m I n_n + k+ \ ik + 



(m 2 + n_k 2 ) 2 

+ (-i) a /. dk + [ „;;;\ 2 e-*+* + - e """ + - . (3.7) 

We must show that M Q , Q ,(x) does not diverge when x + = 0, because in what follows we make 
use of the values of M aa (x — y) at x + = y + to calculate the equal x + -time anticommutation 
relations and to calculate the dynamical operators. 

2 

Since (p + ) 2 = p_ 2 + m 2 , only the term — V among the physical contributions on the first 
line of (|3-7|) gives rise to a linear infrared divergence. We rewrite the expression as follows 

1 r°° dp_ m 2 _ ip _ x - f°° dq m 2 f°° dp_ 1 — cosp^x 



e -v-« = / _| - l^L- / ^- , (3-8) 

4 J-oo p + p_ jo <r 2 Jo p + p_ 

where the first, linearly diverging integration is obtained by changing integration variables 
from p_ to q through p_ = ^ppy- 

As for the residual gauge field's contributions: those divergences which can be regularized 
by the PV prescription we regularize by the PV prescription. Then we find that when x + = 0, 
the first term on the second line of ( |3-7| ) is equal to zero: 



oo rn rn^ roc />^»2 







Tl—Tfl f Tfl 

dk+ m 2 + n.k 2= J dk+ kJ^r 2=0 > (3 ' 9) 



where a = -7==- However, we cannot regularize a linear divergence resulting from a double 
pole by the PV regularization. Therefore when x + = 0, the second term on the second line 
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of (p-7|) gives rise to a linear divergence as follows 



00 „ n_n+m , k+ n+a 



dk+- — . ns „ = — / dk 



00 



+ (m 2 + n„£; + 2 ) 2 4 Jo + \(k + - a) 2 (k+ + a) 2 J 



4 \7o /c 2 ./o k 2 J a k 2 J Jo q 2 2 

where the last, linearly diverging term is given by changing integration variables from k to 
q through k = ^rrq. Thus we see that the linear divergence originated from the physical 
contributions is canceled by that resulting from the residual gauge field's contributions. It 
should be noted here that if we had chosen the representation of the residual gauge fields 
so as to realize the PV prescription (the expression given by ( |2-26| )), we would obtain the 



identically vanishing integral /f^ dk + (^2"^^ ^2 = so that the linear divergence resulting 



from f-X 1 is not canceled. 



It turns out that when x + = 0, the terms on the third line of ( |3-7| ) give rise to a finite 
value, log -, K , so that by summing all the finite contributions we obtain 

y/\K 2 -a 2 \ 



M aa {x)\ x+=0 = K {m \x~\) + i(-l) a - l -e(x-) 
m 2 f°° dp^ 1 — cosp-X~ n + k 



TV 



+ -L + (-!)- log 



2 Jo p+ pj 2 J\K 2 -a 2 



log + log(x- + feC-l)- 1 )) + i{-l) a - l \ + F a (\x-\), (3-11) 

where Ko(rrio\x~\) is the modified Bessel function of order 0, which possesses a logarithmic 
term: —(log +7)- We have gathered the remainder of Kq and the other terms to make 

a continuous function F a (\x~\). 

Now that we have finished showing that the linear divergences cancel each other, we 
proceed to study other properties of the fermion field operators. 

(1) The equal x + -time anticommutation relations hold. For example, when x + = y + , we 
find that 

exp|AM* - y)\ - ^ x -_;^^_ 1) „- 1 e f -°'-»-". (3.12, 



so that if we take to be 



then we obtain 



(z (a) )2 = ?** e -F a ( ) t (3 . 13) 

4:71 



KW,i(!/)}L 

exp[F a (\x- -y-\)-F a (0)} / ij-l)^ 1 ij-l) ' 1 

27r(7°7+) aa \x~ - y- + iei-l) - 1 y~ - x~ + ^(-l)"" 1 
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x expHV^'M - 4,->(!/))]exp[-«Vi(4, +> W - 

When {tp^(x), il>p(y)}\ x +=y+ vanishes due to the Klein transformation factors. 

(2) The vector current, J M = me^d^X, is obtained this time by applying the equal 
x + -time point-splitting procedure to eip^^ip. 

(3) The conserved chiral current J^ 5 = me^ v d v X is similarly obtained from e-07 M 7 5 -?/>. 

(4) As a consequence of our choice of representation of the residual gauge fields, we obtain 
the ML form of the x + -ordered gauge field propagator: 

D^{x -y) = (0\{6(x+ - y + )A,(x)A v (y) + 9(y + - x + )A v (y)A^(x)}\0) 

= (2^)2 / d 2 qD^q)e- iq<X - y) (3-15) 

where 

D^{q) = — 2 j-T- -£/«/ + — — r^--n 1 



q 2 — m? + ie \ q~ + iesgn(g + ) (g_ + 2£sgn(g + )) 2 / 

We give a detailed derivation of this in the Appendix. 

(5) Applying the gauge invariant point-splitting procedure to the fermi products in the 
symmetric energy-momentum tensor and subtracting a diverging c-number, which we denote 
by R, we get 

(9 + ~ = — — i?('07~ d-ip — d-vfj^vp) — R{tp^i~ ^ b d-tp — d-tp^'rfijj) 

- l(A + ,B) + = i(3 + A,cTA) + - l(A + ,B) + , (3-17) 

= -— {(d + A) 2 + (d_A) 2 } + -(F_+) 2 , (3-18) 



^ J R(^7"7 5 5_^ - <9_^7"7 5 ^) + ^^(^7 + 7 5 a_^ - d^ + ^ b ip) 
^{(d + A) 2 + (d_A) 2 } + i(F_ + ) 2 , (3-19) 
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+ = ^ 




(3-20) 



We must take special care in the axial gauge formulation when we derive the conserved 
translational generators from the divergence equation 



The problem is that although x~ is a space coordinate, Q ~ does not vanish in the limits 
£~^±oo. This is because the residual gauge fields do not depend on x~ and because A + 
depends explicitly on x~ . Therefore we have to retain lim. r -^-|- 00 (9 M ~. To take this fact 
into account, we integrate the divergence equation d v Q^ = over a closed surface shown 
in Fig. 1, whose bounds T and L are let tend to oo after all calculations are finished. We 
remark that we can use this surface even in the limit 9 — > | — 0, in contrast with one used 




0. 



(3-21) 



in Refs. 8) and 20). 




X 



X 




Fig. 1. 



It is straightforward to obtain 



= 




-i x =L 



. x =—L 



) 



(3-22) 
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where 



0-(x) 



-ix+=T 



r e+(x)\ x+=T -e+(x) 



x =L 



L O-{x)\ x - =L -0-{x) 



(3-23) 



In what follows we refer to operators made of E as physical operators and operators made 
of B and C as residual gauge operators. 

It can be shown that products of physical operators and residual gauge operators vanish, 
so that the physical operator parts and the residual gauge operator parts decouple in (|3-22|) . 
This result reflects the fact that the physical operator parts are conserved by themselves. It 
can be also shown that the residual gauge operators which are not well-defined in the limit 
L—>oo cancel among themselves. After getting the decoupled and well-defined expressions, 
we take the limit L— >oc. This limit enables us to discard the physical operator parts 6 ~(x) 
of G.r( j cts is usually done when one calculates the conserved physical generators in the 
traditional axial gauge formulations. Finally, we take the limit T— >oo. At this stage we 
assume that we can integrate the remaining residual gauge operators in ~(x) by parts in 
the x + direction. This assumption is justified because the residual gauge operator parts are 
conserved by themselves. As a consequence, we obtain 







dx 



0V=OO 
£ + = — OO 



dx A 



B- 



n-d. 



1 U + 



J x~ = — oo 



(3-24) 



where 



roo 




x +=oo 1 r°° 


/ dx~ 




+ - / dx + 






x+=-oo 2 J -oo 


K + = 


n+ 


: d_Ed + E : 1 : { 




n_ 


2n_ 1 


+ = 


: a_ 


Ed + E : . 



B(x) 



Tl- 



m? — ri-d. 



-,B{x) 



: {(d + Ef + {d-ET + 



(3-25) 



(3-26) 



^From these relations we obtain the conserved generators: 



dx 9jj(x) 



dx + B(x) 



m 2 — n_9 + 2 



/OO j POO 

dx~6_ + (x) + - / dx + B(x 
-oo 2 J — oo 



71- 



d+C(x), 



;B(x) 



(3-27) 
(3-28) 



m 2 — n_9 + 2 

In this way, even in the axial gauge formulation, we can obtain the conserved translational 
generators. These consist of the physical operator part integrated over x + = constant and 
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the residual gauge operator part integrated over x~ = constant. By making use of the com- 
mutation relations Q2-20 ) and ( |3-3| ), it is easy to show that the Heisenberg equations for A + 
and i/j^ hold. This shows that we have successfully constructed the extended Hamiltonian 
formulation of the axial gauge Schwinger model. 

Finally, we shall discuss the charges in the axial gauge formulation. The current, being 
gauge invariant, is the same as in the temporal gauge 



2 2 



J- = -d JmE+— 7^dl l B), J + = d^(mE + — T^d^B) = md.S. (3-29) 

m 2 — nJo% m 2 — 

Using the fact that the current has zero divergence and also using the relation 

m 2 „ „ n_d 2 



)2 B = B+ — 2 B, (3-30) 
m 2 — n_o+ m 2 — n_a\. 

we find the charge to be 

/"OO f-OO _ 

Q = - dx + B(x)+ dx-md-E(x). (3-31) 



OO 



In the temporal gauge formulation (and also in the light-cone gauge 9 -* and Landau gauge 18 ) 
solutions) we set the second term equal to zero. The justification for that is that S(x) is 
a massive field and since it has no infrared singularities, the term involving the integral 
over it commutes with the -0 field. The question is somewhat more complicated in the axial 
gauge case since the singularity coming from dz 1 causes the commutator of the second term 
in ( |3-31| ) with fermi field to become nonzero. On the other hand, we cannot keep both 
terms in ( |3-31| ) since in that case the fermi field would carry twice the correct charge. We 
believe that the correct charge operator is obtained from (|3-31| ) by letting the second term 
provide vanishing topological charges. In previous attempts to formulate the Coulomb gauge 
Schwinger model (using a representation space which is entirely physical), the fermion field 
operators were represented as solitons made solely out of the E so that the E tended to 
nonvanishing values in the imits x~— *-±oo. That led, inevitably, to infrared difficulties. In 
our formulation, the fermion field operators are constituted in the same manner as ones in 
the temporal gauge formulation so that it is reasonable for them not to carry any topological 
charges. To attain this we define the operator (cL) -1 as a finite integral: 

1 r L 



e(x- - y-)f(y-)dy- (3-32) 

2 J-L 

and take the limit L— >oc after all relevant calculations are finished. This procedure is 
justified by the fact that the infrared divergences are all eliminated. We therefore take the 
charge operator to be 

POO 

Q = - / dx + B{x). (3-33) 
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Similarly, we take the chiral charge to be 

/oo 
dx + d+C(x). (3-34) 
-oo 

§4. Concluding remarks 

In this paper we have shown that ^"-independent residual gauge fields can be introduced 
as regulator fields in the pure space-like axial gauge formulation of the Schwinger model 
by extending the Hamiltonian formalism d la McCartor and Robertson. Because we do not 
have consistent pure space-like axial gauge quantization conditions, we have extrapolated the 
solution in the temporal gauge formulation and checked that it is also a solution in the axial 
gauge formulation. One conclusion to be drawn from the work involves the representation 
of the residual gauge fields: We have found that in the axial gauge we must choose the 
representation which specifies the ML prescription in order for the linear infrared divergences 
resulting from the operator (cL)~ 2 to be canceled. For this cancellation to work it is necessary 
for us to change integration variables to ones proportional to n + = sin 29. It follows from this 
that a pure space-like axial gauge formulation in the ordinary coordinates (this is the case of 
9 — | which gives us the Coulomb gauge) has to be defined as the limit 9— »| — 0. Regulating 
the solution in this way involves information off the initial value surface (x° = 0), but in 
any event the solution includes auxiliary fields and the representation space is of indefinite 
metric. 

In contrast, the case of 9 — | (the light-cone axial gauge formulation) can be defined by 
simply setting 9 = | (but the theory still has to be regulated by splitting the fermi products 
off the initial value surf ace). 9 ' 16 ) This is because, by virtue of n 2 = n_ = cos 29 = 0, 
we do not have the linear divergences resulting from (<9_)~ 2 and from the square of the 
inverse of the hyperbolic Laplace operator except for the contact term in the most singular 
component of the gauge field propagator. We have also found that in order for the equal 
x + -time anticommutation relations to be satisfied in the axial gauges, we have to alter the 
representation of the S from the temporal one to the axial one and make a modification to 
the overall normalization of the fermion field operators. 

In contrast to the axial gauges, it seems that in the temporal gauges we can use either 
the representation of the residual gauge fields which specifies the ML prescription or the one 
which specifies the PV prescription. If we use the latter, then both components of the fermi 
field have vanishing vacuum expectation value. 

Now that we have the operator solution in the axial gauge Schwinger model, we can 
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calculate commutation relations of the A^. It turns out that 

[A^x), A v (y)} = i{-g^A(x-y;m 2 ) + (n^d u + n u d IM )d^E(x-y)~n 2 d^d„E(x-y)}, (4-1) 
where A(x; m 2 ) is the commutator function of the free field of mass m and 

E(x) = ^A(x; m 2 ) - - % ~ 5{x + ) + ^±\— 8{x + ). (4-2) 



It remains to be shown whether we can use (44) to obtain consistent pure space-like axial 



gauge quantization conditions which are needed to quantize interacting pure space-like axial 
gauge fields. We leave these tasks for subsequent studies. 



Appendix A 

derivation o/(3-16) 



In this appendix we give a detailed derivation of the x + -ordered gauge field propagator 
D^q) = I d 2 x(0\{6(x + )A ll (x)A u (0) + e{-x + )A v {G)Ap{x)}\Q)e* x . (A-l) 



It is straightforward to show that the contributions from the E and from the residual gauge 
fields are given, respectively, by 

nP / x i m 2 i ( 2n + q + 2 q + 2 \ i 

D p ++ (q) = — — — 2 = ~2 r~— ~9++ + n 2 > A " 2 

q l — m l + is q_ q l — m l + is \ g_ q_ I q_ 



/oo 
dk 




m 2 + n_k 2 \k + — q + — is k + + q + — ie / 



2n + k + m 2 ( % i 
+ Kq-J t-tu , -™ T : ~ + 



(A-3) 



\n 2 k\ + m 2 ) 2 \k + — q + — ie k + + q + — is j 

Here, we have made use of the fact that q 2 = n^q 2 + 2n + q + q- — n^q 2 and that n 2 = n_. 
Note that the explicit x~ dependence gives rise to the factor 8' (?-). Note also that there is no 
on-mass-shell condition for the residual gauge fields so that there remains a k + - integration. 
As a consequence, there arise singularities resulting from the inverse of the hyperbolic Laplace 
operator. It turns out that when we regularize the singularities by PV regularization, the 
integral on the first line of (|A-3|) is well-defined. In fact we can rewrite its integrand as a 
sum of simple poles: 

n^k 2 ( % i \ ri-q 2 ( i 



m 2 + n_k 2 \k + — q + — ie k + + q + — ie J m 2 + n_g + 2 \k + — q + — ie 
i \ n^a ( i i \ 



k + + q + — ie I m 2 + n_q 2 \k + — a k + + a 



(A-4) 
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where a 



—f^=. Direct calculation then gives 

y/-n- to 

r°° I i. 

dk. 



o 



dk A 



ie k + + q + — is / 
( I 1 N 



-7rsgn(g+), 



0, 



(A-5) 



(A-6) 



o \k + — a k + + a / 
where sgn(g + ) is obtained because the k + - integration region is limited to be (0, oo). On the 
other hand, the integral on the second line of ( |A-3| ) yields a linear divergence. We can see 
that by rewriting the integrand as a sum of simple and double poles: 



(m 2 + ri-k , 2 ) 2 \ k-\ 



is 



k + + g H 



IE 



2n + q + m 2 
(m 2 + n_g + 2 ) 2 



x 



q + — IE 

2 

m 



k+ + q^ 

i 



IE 



+ 



ri-q 
+ n+a — 2 



m 



2\2 



k + + a 



n_ m 2 + ri-q 2 \ (k + — a) 2 (k + + a) 



(A-7) 



We can evaluate the integrations of the first and second terms on the right hand with the 
help of ( | A • 5| ) and (|A-6|) . However, we cannot regularize a linear divergence resulting from 
the double pole by the PV prescription. We show below that this linear divergence cancels 
one resulting from the factor {q 2 )~ x of D p ++ in ( A-2 ). For later convenience we rewrite the 
linearly diverging integration in the form 



dk A 



n + 



(k, 



(k + + a) 



dk A 



n + 



/oo n _ 
(A- 
-oo (/_ 



where we have changed the integration variable from k + to g_ 

(|A5| ), flAp and dED into O y ields 



(A; + — a). Substituting 



q 2 — m 2 + ie 



m 2 g + 2 7rsgn(g + )5'(g_) — 2in + q + 5(q_)7rsgn(q + ) 
1 \ 



— m S(q 



dq_ 



where we have made use of the identity 
1 „„ , 1 



q 2 — m 2 + is 
Thus, for the sum of ( 

D ++ (q) 



S'(q- 



2n+g 4 



m 2 + n^q 2 (m 2 + n_g + 2 ) 2 



(A9) 



(A-10) 



and (|A-9|) we obtain 
/ ( 2n + g + 



9++ + 



q z — m' z + ie \ " 1 1 g_ + iesgn(g + ) 

y_ j -co y_ / y_ J 



l q+ nsgn(q + )5'(q_] 



oo "I 

d<2_— 2), 

qj 



(All; 
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where we have made use of the identity: 

2 ~ m * . = %-) f i + ^4tA ■ (A- 12) 

g 2 — m 2 + y q 2 — m 2 + is J 

Now we show that the term — 5(q_)f°^ 00 dq_-p I does not give rise to any divergences 
when we calculate D ++ (x) by substituting ( |A-11|) into ( |A-1|) . To show this we consider the 
integration: 

c 1 ( \ poo 1 \ _ 

d( l-- T~~ - - %~ / <*-—) e ~ iq ~ X 

-oo q 2 — m 2 + is \qt J — oo gl / 

2 , ■ x + — o-: 2/ — ■ ( A ' 13 ) 

qj\q — m 2 + is) m 2 + n_q + J -oo g£ 



-OO 



We can rewrite the second line further as 

e -iq-x~ _ -y y ^oo y 

— 3 — • / '/'/ — ( —T- ~ ■ 



-oo q 2 q 2 — m 2 + is J -oo qj 2 \m 2 + n_q+ q 2 — m 2 + is J 

, / e -iq-x-_ 1 j 2n+g+ + n_g_ \ 

"<?- 5 5 5 ~ ~ 1 7 — 5 9V73 5 3 • (A-14J 

oo y qj q — m 4~ is q_{m 2 + n-q + ){q 2 — m 2 + is) J 

We see that the last integrals diverge at most logarithmically; but logarithmic divergences 
can be regularized by the PV prescription so there arise no divergences in ( |AT4 ). This 
verifies that the following identity holds: 



1 f°° 1 

— + msgn(q + )6'(q^) - S(qJ) / dq_ — 

Q_ J — oo g_ 



Pf4r + msgn(g + )5'(g_) = 1 (A-15) 

qt [q- + zesgn(g + )) z 



where Pf denotes Hadamard's finite part. It should be also noted that 



27r7 -co \ 7 -oo qt qt) nJ o gl 2 

where the Fourier transform of the last term is — 4f? — + ? — ^n)- Thus we obtain 



= -5 T—Z ~9++ + 



2n + g + n q 



2„ 2 
+ 



q 2 — m 2 + is \ g_ + zesgn(g + ) (g_ + ie:sgn(g + )) 2 y 

This completes showing that, due to the residual gauge fields, the linear divergence resulting 
from the d+ l is eliminated from the x + - ordered gauge field propagator. 
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